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, Based on ideas introduced in a previous preprint :ond-mat/9701206 we propose an exactly 

0^ ■ solvable model of bosons interacting amongst themselves via a Van-der Waal-like repulsive 

interaction, and compute both the filling fraction and the dynamical structure factor ex- 
actly. The novelty of this approach involves introducing, analogous to Fermi sea (or surface) 
displacements, Bose fields that in this case, correspond to fluctuations of the Bose conden- 
sate. The exact dynamical structure factor has a coherent part that corresponds to the 
Bogoliubov spectrum and an incoherent part that is a result of correlations. 
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\ I. INTRODUCTION 

The microscopic thoery of He4 is by now a well established and extremely succesful theory. This theory is 
the culmination of the efforts of a large number of distinguished physicists starting from Bogoliubovfl whose 
theory of a weakly non-ideal bose gas was the inspiration for all the theoretical work that followed in this 

3 : 

field. However, the modern theory of superfluidity as it is understood today was poineered by Penrose and 
OnsagerH who introduced the concept of off-diagonal long range order. First attempts at guessing the form 



of the excitation spectrum was made by Bijl and FeynmanEI. Landauu surmised many of the properties of 
the excitation spectrum including the phonon and the roton parts before it was measured experimaentally. 

•rH . 

^ . The microscopic theory for obtaining the static structure factor employs correlated basis functions that was 

pioneered by Bijl, Dingje and Jastrowi. An elaborate application of this technique to quantum liquids in 
general may be found in the book by Feenbergi. Also the work of Kebukawa, Sunakawa and YamasakiB in 
computing the dynamical properties including the excitation spectrum, backflow e.t.c. has to be mentioned. 

In the light of these historical remarks, it seems that yet another theory of this system is redundant. 
However, there is always room for a fresh perspective on well understood problems. The hope is that the 
approach that we are now going to describe along with its fermionic counterpart (described elsewhere! ) 
will produce a model of matter that has these desirable qualities, namely, it is exactly solvable, captures 
the important physics and contains no phcnomenological parameters other than the electric charge and 
mass. If thse types of models were to prove succesful, an outcome that is tantalisingly plausible, it will be a 
vindication of the soundness of the reductionist philosophy. 



1 



II. BOSONS AT ZERO TEMPERATURE 



In this article we do what we did earlier ( |cond-mat/9701206| ) except that here we are bosonizing the 
bosons. Let 6 k and 6^ be the Bose fields in question. Analogous to Fermi sea displacements we introduce 
Bose condensate displacements ( <i k (q)) which are also bose fields obeying canonical commutationrelations 
identical to its parents 6 k and b^. 

& k+q/2 & k-q/2 = N <5 k , <5q,0 + (VN) [5 k +q/2,0 d k(-q) + <5 k _ q/2 , 4(<l)] 
+ Tl ( k ' q ' C l 1 ) d k+q/2-qi/2(qi) d k-q 1 /2(-q + Hi) 

qi 

+ H T 2( k > <J> qi)4-q/2+q 1 /2(qi) d k+q 1 /2(-q + Ql) (1) 

qi 

The above relation is meant to be an operator identity. In other words, all dynamical moments of the 
product fr k+q/ / 2 ^k-q/2 ar e equal when evaluated either in the original bose language or in the condensate- 
displacement language. This is true provided we identify the vacuum of the e? k (q) with the bose condensate 
(the ground state of the noninteracting system). 

d k (q)\BC)=0 (2) 

Furthermore, the chemical potential of the condesate-displacement bosons is taken to be zero. In other words 
the number of condensate displacement bosons is not consereved. 

To make all the dynamical moments of & k+q / 2 ^k-q/2 come out right and the commutation rules amongst 
them also come out right provided we choose T\ and T 2 such that, 

(VA0 2 (5 k+q/2 , <V'-q" /2,o Yl T i( k ''^'' c li)(^k(-q)rf k ' +q ' /2 _ qi /2( t li)rf k '-q 1 /2(- < l' + (O) 

qi 

+ (ViV) 2 5 k+q/2 , V'-qV2,o5] T 2(kW (q")) 



qi 



N <Wq/2,0 ( Wq/ 2 , k "- q "/ 2 ( 5 k -q/ 2 . k ' + q7 2 V-q72. k "+q'72 (^) 

This means, 

T 1 (-q-q72,q',-q) = l; q^O; qV0 (4) 



T 2 (-q/2,q,-q)=0; q ? 0; q ^0 (5) 
in order for the kinetic energy operator to have the form, 
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K = ek d (i/2)k( k ) d (i/2)k(k) (6) 

k 

we must have, 

ek+q^THk + qi /2, 0, qi ) + e k _ qi/2 T 2 (k - qi/2, 0, qi ) = 5 k . qi/2 e qi (7) 

In order for, 

^6 k 6 k = 7V 

k 

we must have, 

Ti(k + qi/2,0,qi)+T 2 (k-qi/2,0,qi) = (8) 

In order for the commutation rules amongst the ^ k+q / 2 ^k-q/2 to come out right, we must have in addition 
to all the above relations, 

T 2 (q/2,q,-q')=0 (9) 

Ta(-q/2,q,q + q) =0 (10) 

Any choice of T\ and T 2 consistent with the above relations should suffice. From the above relations (Eq.( ^) 
and Eq.( ^)), we have quite unambiguously, 

T 1 (k + q 1 /2,0,q 1 )=5 k>qi/2 (11) 

T 2 (k- qi /2,0, qi ) = -5 k>qi/2 (12) 

T 1 (-q-q72,q,-q) = l (13) 

and, 

Ti(k, q, qx) = 0; otherwise (14) 

T 2 (k, q, qi) = 0; otherwise (15) 
This means that we may rewrite the formula for ^ k+q / 2 ^k-q/2 as follows, 

& k+q/2 & k-q/2 = NS^fiSqfi + (VN)(1 - 5 k ,0^q,o) ^k+q/2,0^k(-q) + *k-q/3,0<4(q)] 

+ d (l/2)k+q/2( k + q/ 2 ) d (l/2)k~q/ 2 (k - q/2) 



- <5k,o<5 q ,o £ d qi/2( t * 1 ) d qi/ 2 ( < l 1 ) (16) 
qi 

Consider an interaction of the type, 

H i = (y)E w i£[ <5k +fl/2.o d k(-q) + <5 k _ q / 2 ,o 4( c l)][ <5 k'-q/2,o d k' (q) + V+q/2,0 4'(-q)] ( 17 ) 
q#o k,k' 

or, 

^ = (f ) E %t d -q/2(-q) + < /2 (q)]K/ 2 (q) + dL q/2 (-q)] (18) 

q#0 

and the free case is given by, 

H ° = S e q d (l/2)q( q ) d (V2)q(q) (I 9 ) 

q 

The full hamiltonian may be diagonalised as follows, 

# = £^/q/q (20) 



q 



and, 



A ^ ( " q+ 2 V^ q )^q/2(q) + + ^ + L q/2 (- q) (21) 

/lq ^ ( ^y^ )^q/2(q) + ( ^+^+^q ) ^t q/2( _ q) (22) 

d q/2 (q) = (^±|3±Ma )i/q _ ( -^ + ^ + ^ )*/t q (23) 

Z W q Z CJq 

dt _ q/2 (-q) - ( ^ + 2 V /W W « - C" q+ 2 £q + ^ q )^ q ^ 

^' Z (^ q Z CJ q 



From this we may deduce, 



^ = \/ £ q + 2 /°0 W q £ q ( 25 ) 



( d ; i/2)q (q) d (i/2)q(q)) - ^^X^ (26) 



It is worth while to note that e q + poVq > w q . From this it is possible to write down the filling fraction. 
FILLING FRACTION 

f = N /N = 1 - (l/N) E( d (l/2)q(q) d d/2)q(q)) (27) 
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or, 



f = N /N = l-{l/2^p ) j~ dqqi( ^+Js±^S. } m 



2 UJq 



First define, 

A q = (^±ifL±^)l (29) 

2 Wq 

= + + (30) 

The density operator is, 

p q (t) = ^[d_ (1/2)q (-q)(t)+d| 1/2 ^ (31) 

k 

d-(i/2)q(-q)(*) = ^/-qe-* - B q / q e J (32) 

d (l/ 2)q (q)(*) = ^q/q^ - SqZ-qe"* <*«•* (33) 

d(i/2)k-q/ a (k - q/2)(t) = A k _ q/2 / k _ q/2e - J <■"<-«/>' - B k - q / 2 /! k+q/2 e 4 Wk -/ 2 * (34) 
d (i/2)k +q /2( k + °l/2)W = A k+q/2 /t +q/2 e' - S k+q/2 /_ k _ q/2 e- i <*W (35) 

Now define, 

S>(qt) = ( Pq (t) P _ q (0)) - A r (M-(l/2)q(-q)(t)+4/2)q(q)(*)][ rf (l/2)q(q)(0) +rf t _ (1/2)q (-q)(0)]) 

EKi/2)k +q /2( k + q/2)(t)d(i/2)k-q/a(k - q/2)(t)df 1/2)k ,_ q/a (k' - q/2)(0)d (1/2)k / +q/2 (k' + q/2)(0)) 

k,k' 

= N(— ) ea:p(-i w q i) 
+ £<iW2/-k-q/2e-^/ 2 M k _ q/2 / k ^^ 

k,k' 

+ E^k + q/2/-k-q/ 2 e- i ^+'«/ 2 *S k _ q/2 /t k+q/2e ^-< 1 / 2 *B k ,_ q/2 /_ k , +q/2 B k , +q/2 /^ 
k,k' 

S*>(q,i) = 7V(-^-) «cp(-i w q i) 

+ ^2eXp(-l (^k+q/2 + W k-q/2)0[ B k+q/2^k-q/2 + £k+ q /2-B-k+q/2-4k- q /2-4k+q/2] 



and, 



S < (<1, t) = N (—) exp{i ujqt) + V exp(i (w k+q/2 + w k _ q/2 )i) 
w o k 

[ B k-q/2 A k+q/2 + -Bk-q/2-B_ k -q/2^k+q/2^k-q/2] 

S>(q,t) = (p q (t)p- q (0)) 



S<(q,t) = (p- q (0)p q (t)) 

From the above equations, it is easy to see that there is a coherent part corresponding to the Bogoliubov 
spectrum and an incoherent part which is due to correlations and is responsible (hopefully) for the roton 
minimum. 
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